Where do perturbative and non-perturbative QCD meet? 
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We computed the static potential and Wilson loops to 0(a 2 ) in perturbation theory for different lattice quark 
and gluon actions. In general, we find short distance lattice data to be well described by "boosted perturbation 
theory". For Wilson- type fermions at present-day quark masses and lattice spacings agreement within 10 % 
between measured "/3-shifts" and those predicted by perturbation theory is found. We comment on prospects for 
a determination of the real world QCD running coupling. 



1. INTRODUCTION 

Perturbative results on physical quantities are 
useful in lattice simulations in several ways: 

• to qualitatively understand "improvement" and 
"lattice artefacts" , 

• to numerically validate that the continuum limit 
exists and is reached as (3 — > oo, 

• for predicting the "/3-shift" , resulting from mas- 
sive sea quarks, 

• to estimate quantities that are hard to obtain 
otherwise, e.g. the static quark self energy, 

• to calculate a s from low energy QCD phe- 
nomenology. 

In view of this we calculated Wilson loops 
and the static potential with massive Wil- 
son, Sheikholeslami-Wohlert (SW) and Kogut- 
Susskind (KS) fermions to 0(a 2 ) f). We de- 
fine the potential, aVS nt (Ra) = VL,i(R)ctL + 
v Lt2 (R)a 2 L + where Vj nt (Ra) = V(Ra) - 
U(oo). The static quark self energy, U(oo) = 
2 (5?7j s tat oc aa^ 1 + • ■ •, vanishes in dimensional 
regularisation but diverges on the lattice as the 
continuum limit is approached. For SU(3) we 
find, aSm stat = 2.1172743 a L + [11.143(3) + 
n s Y s \a\ with Y f = -0.36846(6) for massless KS 
quarks and Y f = -0.42333(6) + 0.0516(2)c sw - 
0.5870(2)c| w for Wilson-SW quarks. While our 
value of the Wilson-SW Yf agrees with Ref. @ we 
find a 4 (j discrepancy in the gluonic contribution. 
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Figure 1. Wilson quark contribution to vl^- 

2. "AKi" AND THE "/3-SHIFT" 

The MS scheme is related to the lattice scheme 
via ajgg(a _1 ) = + b\o? + ■ ■ •, with the con- 
version factor &i = — ir/(2N) + k\N + Ki(ma)rif, 
where the numerical constant k\ is known for a 
variety of gluonic actions and Ki(0) is known 
for Wilson, SW and KS quarks with Wilson 
glue. Further couplings can be defined, e.g. 
from the potential in position space, cxr((J.) = 
-rV int (r)/C F , fi = r _1 . While the ~M~S scheme 
is "mass-independent" , the /3-function coefficients 
Pi of the above "i?" scheme depend on mj \i and 
hence on the quark mass; universality is lost and 
the (3i are specific for each dimensionful observ- 
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Figure 2. Matching with the MS scheme. 
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Figure 3. Comparison with the measured /3-shift. 

able that is studied. The same holds true for the 
conversion factor to the MS scheme. 

In contrast, the lattice scheme is mass- 
independent in the continuum limit. However, 
it is often worthwhile to study quantities that 
are not defined in this limit like small Wilson 
loops ||, and some effective field theory ap- 
proaches require a finite lattice cut-off too. At 
fixed physical r = Ra the limit R — > oo corre- 
sponds to the continuum limit a — > 0, in which 
rotational symmetry is restored and lattice and 
continuum perturbative predictions agree. One 



would also expect rotational symmetry to be re- 
stored at finite a for distances r ^S> a. These 
two cases become distinguishable in perturbation 
theory once an external scale m is introduced and 
hence the situation r < m differs from r > m. 

At finite a and m > the lattice scheme be- 
comes mass-dependent too, as indicated by the 
function K\(md) — Ki(0) + AK~i(ma) within b\ 
above. We denote the fermionic contribution to 
vl,2 as rifv^ 2 and V[ 2 — w{ 2 /(47r) 2 . The lat- 
ter, multiplied by —R, is displayed in Fig. |f| for 
the example of Wilson fermions. In the massless 
case we find the expected logarithmic running, 
proportional to the fermionic contribution to /3o 
(straight line). The offset at small R is related 
to K\ (ma) , which we determine by matching the 
lattice potential at R ^> 1 to the known M S re- 
sult. Aifi is only unique in the limit m = 
(where it vanishes). For massive quarks univer- 
sality is lost and AK~i will depend on the observ- 
able that is matched. In the case of the position 
space potential (and force) we obtain the asymp- 
totic behaviour, AK 1 (ma) = 0.0011(3) - A^(0) - 
ln(ma)/(37r) for ma —> oo. The logarithmic term 
is universal and guarantees the massive fermions 
to decouple from the (3 function. 

We display the size of this correction, rela- 
tive to Ai(0) < in Fig. |. Note that in the 
limit ma — > oo the ratio — AK\/K\ diverges to- 
wards negative values. For the two O(a) im- 
proved fermionic actions the mass dependence is 
minimal while this effect can be very significant 
for Wilson fermions and has to be taken into ac- 
count in any calculation that uses the potential 
as an intermediate scheme. 

As can be seen from Fig. [lj at r <C m _1 , i.e. 
R <C (ma) -1 , the effective Coulomb coupling is 
screened with the same logarithmic slope as in 
the massless case. At r ~S> m~ x the heavy quarks 
decouple and do not contribute to the running 
anymore but just to the overall normalisation. 
It is here that the un-quenched potential can be 
matched to the quenched one by adjusting the 
coupling constant (3 = 3/(27ra L ): = 0^ + 

A/3, A/3 = n s [ln(Dma) + 37rAAi(ma)] /(2tt 2 ) + 
0((3~ r ). We find the numerical values, D = 
0.448(2), 0.0238(1) and 0.726(2), for Wilson, SW 
and KS fermions, respectively. As ma — > oo the 
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AKi term guarantees that A/3 — > while A/3 di- 
verges as ma — > 0: in the light quark limit the 
quenched and un-quenched theories decouple. 

In Fig. H we compare our perturbative predic- 
tion to numerical data with n/ = 2 Wilson [Q, 
SW !| and KS |6) fermions, obtained at lat- 
tice spacings 4 < ro/a < 6.5 and quark masses 
0.1 > ma > 0.01. The Wilson and SW results 
fall onto a universal curve that differs by less than 
10 % from the prediction while the KS results de- 
viate much more and some /3 dependence is evi- 
dent. Whether this is due to a slower convergence 
of the perturbative series or due to n/ not being a 
multiple of four is an open question. The qualita- 
tive agreement between prediction and simulation 
for Wilson- type quarks indicates that, at least at 
present masses, physics at hadronic scales is not 
strongly affected by quark loops, which is con- 
sistent with the phenomenological success of the 
quenched approximation. 

3. THE QCD RUNNING COUPLING 

We determine the running coupling 
ctjjg(3A01 a -1 ) from the average plaquette, fol- 
lowing the method detailed in Ref. [3]. The result 
is then converted into units of r |Q- 6 | and numer- 
ically evolved to the scale fi = 10 (Fig. |^) via 
the four-loop (3 function. While the mass depen- 
dence of the plaquette is rather weak, the term 
AKi(ma) has a statistically significant impact 
on the Wilson results at masses ma > 0.02. 

The bottom-left n/ = pentagon corre- 
sponds to the ALPHA collaboration result, 
a-g-g(n) = 0.157(4), the other ones have been 
obtained from the average plaquette within a /3- 
window that covers the range of lattice spacings, 
spanned by the dynamical simulations. We note 
that the n f — values obtained by use of the two 
different methods agree. All dynamical results 
can be fitted to polynomials in mro (solid curves). 
The SW results || have been obtained at fixed 
r while the two Wilson || and three KS || 
data sets have been produced at fixed j3. We 
observe a significant lattice spacing dependence, 
in contrast to the nj — case. Extrapolating 
to the limits mro ~ 0.01 and a/ro — » 0, the 
latter linearly for Wilson and quadratically for 
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Figure 4. QCD coupling for rif = 2. 

KS fermions, we obtain the consistent estimate, 
a^L(n) - a^L(^i) = 0.030(6). Simulations at 
additional lattice spacings are mandatory to con- 
trol the extrapolation to the continuum limit. 
Furthermore, numerical data at n/ ^ 0,2 are 
essential, before a^-^{mz) can be predicted with 
confidence. Once such results are available an 
error of 2 % appears to be realistic. 
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